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We compute the electromagnetic form factors of the nucleon in the Poincare-covariant Faddeev 
framework based on the Dyson-Schwinger equations of QCD. The general expression for a baryon's 
electromagnetic current in terms of three interacting dressed quarks is derived. Upon employing a 
rainbow-ladder gluon-exchange kernel for the quark-quark interaction, the nucleon's Faddeev am- 
plitude and electromagnetic form factors are computed without any further truncations or model 
assumptions. The form factor results show clear evidence of missing pion-cloud effects below a 
photon momentum transfer of ~ 2 GeV 2 and in the chiral region whereas they agree well with 
experimental data at higher photon momenta. Thus, the approach reflects the properties of the 
nucleon's quark core. 
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I. INTRODUCTION 

The measurement of the nucleon's elastic electromag- 
netic form factors has been subject to major experimental 
progress during the past years; see Refs. [1-3] for recent 
reviews. Precise polarization experiments at Jefferson 
Lab, MIT-Bates and MAMI have charted the evolution 
of the electromagnetic form factors over a wide range of 
photon momenta, and the forthcoming 12-GeV upgrade 
at JLAB will allow for an extension of these measure- 
ments even farther beyond. In connection to this, the N* 
program at JLAB is expected to provide substantial in- 
sight in the internal structure of nucleon resonances [4, 5] . 

As electromagnetic form factors encode the nucleon's 
spatial charge and magnetization distribution, they pro- 
vide a direct connection to the underlying dynamics of 
quarks and gluons which are the fundamental degrees of 
freedom in Quantum Chromodynamics (QCD). The new 
data have raised many new questions and stimulated a 
parallel response from theory. Among the aspects that 
have been intensely discussed are the question of two- 
photon contributions to elastic cross sections [6-8]; the 
importance of quark orbital angular momentum in the 
nucleon's structure at high momentum [9-12]; or the role 
of pseudoscalar meson clouds and the interpretation of 
form factors through charge and magnetization densities 
in different frames [13-17]. At the same time, advances 
in studying hard exclusive processes have provided in- 
sight into the spatial and spin structure of the nucleon 
through the framework of generalized parton distribu- 
tions [18-21]. 

Approaching these issues from a fundamental theoret- 
ical perspective aims at a quantitative understanding of 
the hadrons' substructure in terms of quarks and gluons 
in QCD. A comprehensive theoretical framework which 
has seen notable progress in recent years is given by the 
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Dyson-Schwinger equations (DSEs) of QCD. They form 
a set of infinitely many coupled integral equations for 
QCD's Green functions that grants access to nonper- 
turbative phenomena such as dynamical chiral symme- 
try breaking, confinement, and the formation of bound 
states; see [22-24] for reviews. 

The bridge between Dyson-Schwinger equations and 
hadron properties is provided by covariant bound-state 
equations, see [25-27] and references therein. In the case 
of baryons, it is the covariant Faddeev equation [28- 
31] that describes a baryon's amplitude through all in- 
teractions that can take place between its dressed va- 
lence quarks, and thereby constitutes the direct ana- 
logue to the Bethe-Salpeter equation for mesons [32]. 
The Faddeev equation was recently solved for the nu- 
cleon mass [33, 34] with an interaction kernel that traces 
the nucleon's binding to its quark-quark correlations and 
thereby relates the properties of baryons to those of 
mesons and the underlying substructure in QCD. First 
steps toward the A— baryon in that setup were reported 
in Ref. [35], and the heavy-quark limit of the Faddeev 
equation was recently studied in Coulomb gauge [36]. 

In order to compute electromagnetic form factors one 
must additionally specify how the photon couples to the 
dressed quarks in a baryon. Due to the complexity 
of the three-quark problem, baryon form-factor calcula- 
tions within the Dyson-Schwinger framework have been 
performed in the covariant quark-diquark model [37-47] 
which traces the nucleon's binding and its interaction 
with the photon to scalar and axial-vector diquark cor- 
relations in its interior. In Refs. [48, 49] the quark- 
diquark model was extended insofar as the dynamics 
of the quarks and diquarks were determined from their 
quark and gluon constituents, and corresponding results 
for nucleon and A— baryon form factors were reported 
in Refs. [27, 50, 51]. These studies provide valuable in- 
sight in the correlations that take place amongst the nu- 
cleon's constituents when probed by a photon, and sev- 
eral characteristic features revealed in experiments are 
readily reproduced. Nevertheless, the predictive power of 
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FIG. 1: (Color online) The covariant three-quark equation for a baryon amplitude, Eq. (1). 



the quark-diquark approach is limited by the sensitivity 
to the diquark parameters at larger photon momentum. 

In the present work we follow Ref. [33] and solve the 
covariant Faddeev equation directly without introducing 
explicit diquark degrees of freedom. The nucleon's bind- 
ing is described by a rainbow-ladder (RL) interaction 
which amounts to an iterated dressed gluon exchange be- 
tween any two quarks. The single parameter of the ap- 
proach is a scale which is fixed to meson properties. We 
construct a consistent electromagnetic current operator 
and compute the nucleon's electromagnetic form factors 
in this setup. A rainbow-ladder truncation excludes the 
presence of meson-cloud corrections which have signifi- 
cant impact on the chiral and low-momentum structure 
of hadrons, hence the approach reflects the properties of 
a hadronic quark core. 

The manuscript is organized as follows: in Section II 
we recall the Poincare-covariant Faddeev equation, col- 
lect its ingredients, and discuss the nucleon amplitude 
and the current-mass evolution of its mass. In Section III 
we construct the nucleon's electromagnetic current op- 
erator in the three-quark framework. Results for the 
nucleon's electromagnetic form factors are presented in 
Section IV. Technical details of the calculation are rele- 
gated to the appendices: conventions and formulas are 
collected in App. A; the structure of the nucleon ampli- 
tude, its permutation-group properties, basis decomposi- 
tion and partial- wave analysis are analyized in App. B; 
App. C presents an updated solution strategy for the Fad- 
deev equation; and the nucleon's electromagnetic current 
operator and its ingredients in the three-quark frame- 
work are described in App. D. Throughout this paper we 
work in Euclidean momentum space and use the isospin- 
symmetric limit m u = ma- 



il. COVARIANT FADDEEV-EQUATION 

Baryons in QCD appear as poles in the three-quark 
scattering matrix. The residue at a particular pole de- 
fines the baryon's bound-state amplitude. It satisfies a 
covariant homogeneous integral equation which is illus- 
trated in Fig. 1 and reads 



* = K * . 



K = K h 



(1) 



a=l 



where 'J is the bound-state amplitude defined on the 
baryon mass shell. The three-body kernel K comprises 



a three-quark irreducible contribution and the sum of 
permuted two-quark kernels whose quark-antiquark ana- 
logues appear in a meson BSE, and the subscript a de- 
notes the respective spectator quark. 

It has been a longstanding suggestion that the strong 
attraction between two quarks in a SU(3)c antitriplet 
configuration within a baryon is a key feature for a better 
understanding of hadron properties [52, 53]. This has 
been the guiding idea in the quark-diquark approach, and 
in our context it motivates the omission of the three-body 
irreducible contribution, i.e. the first diagram on the 
right-hand side of Fig. 1, from the full three-quark kernel. 
The resulting covariant Faddeev equation attributes the 
binding mechanism of three quarks in a baryon to its 
quark-quark correlations and reads: 



= J [irg, 77 , ^ 7 , 5 ( P w, (7 ( 1 ),p) + 



(2) 



where K ( a ) denotes the renormalization-group invariant 
product of a qq kernel and two dressed quark propaga- 
tors: 



A' 



(a) 



, = K,. aa »pp" S a n a , (p b - k) Spi>p (p c + k) , (3) 



and {a, 6, c} is an even permutation of {1, 2, 3}. 

The notation in Eqs. (2-3) is understood as follows: 
the Poincare-covariant nucleon amplitude carries three 
spinor indices {a, /3, 7} for the valence quarks and one in- 
dex 6 for the spin-1/2 nucleon. It depends on three quark 
momenta pi, P2, P3 which can be expressed in terms of 
the total nucleon momentum P, where P 2 — —M 2 is 
fixed, and two relative Jacobi momenta p and q. Upon 
choosing symmetric momentum partitioning they are re- 
lated via: 
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The quark propagators S depend on the internal quark 
momenta p^ — k andpj+fc, where k is the loop momentum 
that will later be identified with the gluon momentum, 
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FIG. 2: (Color online) Quark DSE (7) in rainbow-ladder trun- 
cation. 



and the internal relative momenta in Eq. (2) are given 
by: 

=p + k, <?W =q-k/2, 

pW=p-k, qW=q-k/2, (5) 

p^ =p, = q + k. 

The inherent color structure leads to a prefactor 2 /3 for 
the integral in Eq. (2). In the following subsection we 
collect the ingredients of Eqs. (2-3), i.e. the dressed- 
quark propagator S and the quark-quark kernel /C. 

A. Quark propagator and qq kernel 

The fundamental ingredient which appears in the co- 
variant Faddeev equation and provides a link between 
hadron properties and the underlying structure of QCD 
is the dressed quark propagator S(p). It is expressed in 
terms of two scalar functions, the quark wave-function 
renormalization l/A(p 2 ) and the quark mass function 

M(p 2 y. 

S- 1 (p)^A(p 2 ) (i# + M(p 2 )). (6) 

Dynamical chiral symmetry breaking is manifest in a 
non-perturbative enhancement of both dressing functions 
M(p 2 ) and A(p 2 ) at small momenta which indicates the 
dynamical generation of a large constituent-quark mass. 
Such an enhancement emerges in the solution of the 
quark DSE, cf. Fig. 2: 

S ap(P) = Z 2 (¥ + m ) a p+ J fc aa >p>p(p,q) S a >f3>(q) , (7) 

Q 

where Zi is the quark renormalization constant. The 
bare current-quark mass m constitutes an input of the 
equation and can be readily varied from the chiral limit 
up to the heavy-quark regime. The interaction kernel K, 
includes the dressed gluon propagator as well as one bare 
and one dressed quark-gluon vertex. 

While in principle the dressed gluon propagator and 
quark-gluon vertex could be obtained as solutions of the 
infinite coupled tower of QCD's DSEs, practical numer- 
ical studies rely upon truncations, combined with sub- 
stantiated ansatze for Green functions which are not 
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FIG. 3: (Color online) Effective coupling a(k 2 ) of Eq. (9), 
evaluated for A = 0.72 GeV corresponding to the u/d— quark 
mass and in the range r\ £ [1.6, 2.0]. 



solved for explicitly. The pion's nature as the Goldstone 
boson of spontaneous chiral symmetry breaking necessi- 
tates a truncation that preserves the axial-vector Ward- 
Takahashi identity. The latter relates the kernel of the 
quark DSE with that of a meson BSE, ensures a massless 
pion in the chiral limit and leads to a generalized Gell- 
Mann-Oakes-Renner relation [54, 55]. Such a symmetry- 
preserving truncation scheme was described in [56, 57], 
and its lowest order is the rainbow-ladder (RL) trun- 
cation which amounts to an iterated dressed-gluon ex- 
change between quark and antiquark. In rainbow-ladder 
only the vector part ~ 7 M of the quark-gluon vertex is 
retained whose non-perturbative dressing, together with 
that of the gluon propagator, is absorbed into an effec- 
tive coupling a(k 2 ) which is modeled. The kernel /C of 
both quark DSE and Faddeev equation then reads: 

2 47TQ(fc 2 ) 

= ^2 p k 7 QQ ' 7/3/3' , (») 

where T£ v = — k^k v is a transverse projector with 
respect to the gluon momentum k and a hat denotes a 
normalized 4-vector. 

At large gluon momenta the effective coupling a(k 2 ) is 
constrained by perturbative QCD, whereas the details of 
the interaction in the deep infrared are not important for 
hadronic ground states [58]. At small and intermediate 
momenta it must exhibit sufficient strength to allow for 
dynamical chiral symmetry breaking and the dynamical 
generation of a constituent-quark mass scale. We employ 
the frequently used ansatz [59] 

a{k 2 ) = w(^) e"" 2 ^) +auv(fc 2 ), (9) 

where the second term reproduces the logarithmic de- 
crease of QCD's perturbative running coupling and van- 
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ishes at k 2 = 0. The first term supplies the necessary in- 
frared strength and is characterized by two parameters: 
an infrared scale A and a dimensionless width parameter 
■q, cf. Fig. 3. 

In combination with the interaction of Eq. (9) , the RL 
truncation has been extensively used in Dyson-Schwingcr 
studies of hadrons. Upon setting the scale A = 0.72 GeV 
to reproduce the experimental pion decay constant, RL 
provides a reasonable description of pseudoscalar-meson, 
vector-meson, nucleon and A ground-state properties, see 
e.g. [60-63] and references therein. Moreover, these ob- 
servables have turned out to be largely insensitive to the 
shape of the coupling in the infrared [59, 62]; i.e., to a 
variation of the parameter rj around the value rj rs 1.8. 

Progress has also been made for other meson quan- 
tum numbers such as axial- vector and pseudoscalar isos- 
inglet mesons whose properties are subject to substan- 
tial corrections beyond rainbow-ladder [64-70]. However, 
such analyses typically require a significant amplification 
of numerical effort which is not yet feasible for stud- 
ies in the baryon sector. Important attractive contribu- 
tions beyond RL come from a pseudoscalar meson-cloud 
which augments the 'quark core' of dynamically gener- 
ated hadron observables in the chiral regime and van- 
ishes with increasing current-quark mass. Such effects 
are missing in a RL truncation. Thus, the present work 
aims at investigating the electromagnetic form factors of 
the nucleon's quark core. 

We note that through Eq. (9) all parameters of the 
interaction a(k 2 ), and thereby all equations that appear 
in subsequent sections, are fixed by using information 
from pion properties only. 

B. Nucleon amplitude 

Upon having determined the input of the covariant 
Faddeev equation through Eqs. (6-9) one can proceed 
with its solution. First results for the nucleon's mass 
and bound-state amplitude were reported in Ref. [33, 34]. 
While the mass can be reliably determined with rela- 
tively modest numerical accuracy, the form-factor com- 
putation requires a significantly higher resolution of the 
Faddeev amplitude, especially at larger photon momen- 
tum transfer. In the present work we use a solution tech- 
nique that exploits the permutation-group properties of 
the amplitude. This enables us to drastically reduce the 
involved CPU times and solve the Faddeev equation with 
its full momentum dependence. The method is described 
in App. C. 

The structure of the nucleon amplitude, together with 
its basis decomposition and permutation-group proper- 
ties, is discussed in detail in App. B. In the following we 
will highlight some key aspects. The spin-flavor structure 
of the on-shell nucleon amplitude can be expressed as 

2 

* = *-F = ^$ n F„, (10) 

71=1 



where the Dirac amplitude \t and the isospin- 1 ^ flavor 
tensor F of Eq. (B2) transform as doublets under the 
permutation group S 3 , with mixed- antisymmetric entries 
Fi and mixed-symmetric components ^2, F2, respec- 
tively. The structure of Eq. (10) ensures the Pauli princi- 
ple: the nucleon amplitude involving its full spin-flavor- 
color structure must be antisymmetric under quark ex- 
change, hence its spin-flavor part has to be symmetric. 
The spinor parts ^ n involve 64 covariant, orthogonal and 
momentum-dependent Dirac structures X^^, 

V n (p, q,P) = Y^ fn,kiju X k,ij^ , (H) 

kijcj 

which are discussed in detail in App. B2. The ampli- 
tude dressing functions f n ,kijuj depend on the 5 Lorentz- 
invariant momentum variables 

P 2 , Q 2 , z =p^-qf, z x =p-P, z 2 =q-P, (12) 

where a hat denotes a normalized 4- vector and the sub- 
script 'T' a transverse projection with respect to the nu- 
cleon momentum P. The total momentum-squared is 
fixed: P 2 = -M 2 . 

The orthogonal basis elements ~Kk,iju are eigenstates 
of total quark spin and orbital angular momentum in 
the nucleon's rest frame; the corresponding partial-wave 
decomposition is explained in App. B3. The rest-frame 
nucleon amplitude is dominated by s— wave components, 
i.e. by the subset of eight relative-momentum indepen- 
dent basis elements which carry total quark spin s = 1 /% 
and orbital angular momentum I = 0. We denote them 
here by 



s± 


; — X^iii 


= A±7 5 C ® A+ , 


v± 


: = Xi^ii 


= 75 7? A ±75 C ® 7? A+ 


p± 


: = X2,ll± 


= (75 ® 75) s± , 


A± 


:= X2 j 2i± 


= (75 ® 75) v± , 



(13) 



where the 7— matrices 7^ are transverse with respect to 
the nucleon momentum P. The remaining basis elements 
are either p— or d— waves. Table I shows the s— , p— and 
d— wave contributions to the nucleon's canonical normal- 
ization integral (23) at different current-quark masses. 
The s— wave elements contribute roughly 2 /3 to the norm 
and the p— waves the remaining third. The p— wave con- 
tribution decreases, albeit very slowly, with higher quark 
mass which signals a substantial amount of orbital an- 
gular momentum in the nucleon's rest-frame amplitude 
well beyond the strange-quark mass. 

To analyze the s— wave components in the nucleon am- 
plitude in more detail it is instructive to rearrange the 
eight basis elements of Eq. (13) in permutation-group 
multiplcts. This yields the orthonormal doublets 

^3^2V-A_j' 

(14) 

73 V A- + V + V3 P J ' 
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FIG. 4: (Color online) Result for the three dominant s— wave 
contributions in the nucleon's Faddeev amplitude. The plot 
shows the zeroth Legendre and Chebyshev moments (in the 
variables j/i , zo and zi, Z2, respectively) of the dressing func- 
tions Si defined via Eqs. (14) and (16). 



where the upper entries are mixed-antisymmetric with 
respect to the first two Dirac indices and the lower entries 
mixed-symmetric, and two further singlets 



4-, 



1 

71 



A_ +V- V3P 



(15) 



which are fully antisymmetric or symmetric, respectively. 
Here we defined P := (P+ + P_)/2, P := (P+ - P_)/2 
and accordingly for V and V. 

Eqs. (14-15) imply that, without including a depen- 
dence on the relative momenta, only three fully sym- 
metric Dirac-flavor combinations 1 jK < ) • F can arise in the 
s = 1 /2, I = subspace. They appear in combination 
with symmetric singlet dressing functions which are lin- 
ear combinations of those associated with the basis ele- 
ments in Eq. (13) and must depend on symmetric combi- 
nations of the momentum variables in Eq. (12). Denot- 
ing them by a fully symmetric spin/momentum-flavor 
amplitude is then obtained via 



(16) 



where the dots refer to further combinations of Eqs. (14- 
15) with mixed- (anti-)symmetric dressing functions, and 
also to the remaining p— and d— wave components. In 
Eq. (B7) we define momentum variables that transform 
as multiplets under S 3 , namely a symmetric singlet vari- 
able 
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TABLE I: s~, p— and d— wave contributions to the nucleon's 
canonical normalization at three pion masses, expressed as 
fractions of 1. The first column corresponds to the physical 
u/d— quark mass. 



and four dimensionless angular variables y\, j/2, wi, W2 
which form doublets. The dressing functions Si can then 
only depend on the variable x and the symmetric combi- 



nations 



■vh w i 



w%, and yiwi + y 2 w 2 - 



(17) 



The full solution of the Faddeev equation indeed re- 
veals the three singlet dressing functions Sj to contribute 
the bulk to the s— wave fraction in the normalization. 
Their angular dependence is weak, especially in the vari- 
ables z-i and Zo, and a corresponding polynomial expan- 
sion vanishes rapidly. The zeroth angular moments of 
the three Sj are plotted in Fig. (4) as a function of the 
variable \fx. All three dressing functions turn out to be 
large; in particular, s\ and S3 are almost identical in size. 

The resulting current-mass evolution of the nucleon's 
mass is displayed in Fig. 5. The pion mass was ob- 
tained from its pseudoscalar-meson Bethe-Salpeter equa- 
tion with the same rainbow-ladder input. The scale 
A in Eq. (9) was fixed to reproduce the experimental 
pion decay constant. In agreement with previous meson 
and quark-diquark studies, the sensitivity to the infrared 
shape of the effective coupling a(k 2 ) is small; this is in- 
dicated by the band which corresponds to a variation 
77 = 1.8 ± 0.2, cf. Fig. 3. At the physical u/d— quark 
mass, our result Mm = 0.94 GeV is in excellent agree- 
ment with the experimental value, and its current-mass 
evolution compares reasonably well with lattice data at 
higher quark masses. 

In connection with Fig. 5 we reiterate that contri- 
butions from a pseudoscalar-meson cloud are absent 
in a rainbow-ladder truncation; the current approach 
can therefore be viewed to describe a hadronic quark 
core. Such corrections can be estimated from chiral 
effective field theory and would yield a reduction of 
~ 20 — 30% of the nucleon's core mass in the chiral re- 
gion [63]. The proximity between our calculated mass 
and the experimental and lattice values therefore sug- 
gests a non-perturbative cancelation mechanism beyond 
rainbow-ladder. Indeed, such a behavior emerges for 
the p— meson where attractive pion-cloud effects beyond 
RL are essentially saturated by further repulsive con- 
tributions from the quark-gluon vertex and the quark- 
antiquark kernel [67-70]. In addition, the second type of 
corrections dominates in scalar and axial-vector mesons 
which explains why these quantum numbers arc not 
well reproduced in a RL truncation. Given the quali- 
tatively quite similar behavior of the p— meson mass in 
the present framework in comparison with lattice data, as 
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FIG. 5: (Color online) Current-mass evolution of the nu- 
cleon mass obtained from the covariant Faddeev equation. 
The band corresponds to a variation of the width parameter 
rj = 1.8 ± 0.2; the star denotes the experimental value. We 
compare to a selection of lattice data from Refs. [71-77]. 



well as nucleon and A masses in a quark-diquark calcu- 
lation [49, 78], it is conceivable that a similar mechanism 
might be at work here. 

We note that the Faddeev result for the nucleon 
mass in Fig. (5) is in excellent agreement with that ob- 
tained from the quark-diquark calculation (cf. Fig. 2 in 
Ref. [78]), where the qq scattering matrix was reduced to 
a sum of scalar and axial-vector diquark channels which 
were solved from their Bethe-Salpeter equations. This 
highlights the importance of the diquark concept as rel- 
evant degrees of freedom within hadrons, and it specifi- 
cally shows that including only the lightest diquarks can 
already account for the bulk of the nucleon mass. In Sec- 
tion IV it will turn out that this observation also holds 
for the overall behavior of the nucleon's electromagnetic 
form factors. 



III. ELECTROMAGNETIC CURRENT 

A. Construction of the current 

Having numerically computed the nucleon amplitude, 
we proceed with the construction of the nucleon's elec- 
tromagnetic current. Its general structure is given by (cf. 



j*(p, Q) = iki [f^ - A a^cy) a; (18) 

and parameterized by two form factors: the Dirac form 
factor Fi(Q 2 ) and the Pauli form factor F 2 (Q 2 ). If P, 
and Pf denote the initial and final nucleon momenta, the 
exchanged photon momentum is given by Q = Pf — Pi 
and the average total momentum by P = (Pi + Pf)/2. 
The positive-energy projectors A^' 1 = k + (Pj.i) are de- 
fined in Eq. (A10). The Dirac and Pauli form factors can 
be expressed through the Sachs form factors: 



Gi 



F 1 - tF 2 



G 



M 



F\ + F<i , 



(19) 



with r = Q 2 /(AM 2 ). Their static dimensionless values 
are the proton's and neutron's charges G^ n (0) = {1,0} 
and their magnetic dipole moments G^ n (Q) = /i p ' n . 

In order to relate the electromagnetic current of 
Eq. (18) to the underlying description of the nucleon as a 
composite object, one must specify how a photon couples 
to its constituents. A systematic construction principle 
for the nucleon-photon current based on electromagnetic 
gauge invariance is the 'gauging of equations' prescrip- 
tion [79-81] which was previously used to derive the elec- 
tromagnetic current in the quark-diquark model [42, 82]. 
We will briefly sketch the procedure here. Gauging, in- 
dicated by an index fj,, has the properties of a derivative, 
i.e. it is linear and satisfies Leibniz' rule. In the present 
context, fi denotes the coupling to a photon with mo- 
mentum Q: when acting upon an n— point Green func- 
tion, it yields an (n + 1)— point function with an addi- 
tional photon leg. For instance, the quark-photon vertex 
is the gauged inverse quark propagator: = (5 f_1 ) AI , 
and hence = -S (S^ 1 )^ 1 S = ~ST% S. 

To derive a baryon's electromagnetic current in this 
framework, we recall the nonperturbatively resummed 
Dyson series for the three-body scattering matrix T, 
i.e. Dyson's equation, which is the starting point of the 
bound-state approach: 



T = K + KG T 



T- 1 = K- 1 - Gn 



(20) 



where Go is the product of three dressed quark propaga- 
tors and K = K Gq 1 is the three-body kernel of Eq. (1). 
A pole in the scattering matrix defines a bound state on 
its mass shell P 2 — —M 2 . The scattering matrix at the 
pole assumes the form 



T 



P 2 + M 2 



(21) 



where '5 is the baryon's bound-state amplitude and ^ its 
charge conjugate, and we suppressed all indices as well as 
the momentum dependence and dimensionful prefactors. 
At the pole, Dyson's equation reduces to the homoge- 
neous bound-state equation (1) for the amplitude \&: 



K G * = * 



T- 1 ^ = 0. 



(22) 
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FIG. 6: (Color online) General expression for a baryon's electromagnetic current (modulo permutations) in the three-quark 
framework. 



Evaluating the relation T = -T (T' 1 )' T at the bound- 
state pole, where ' denotes the derivative d/dP 2 , yields 
in combination with Eq. (21) the on-shell canonical nor- 
malization condition: 



yields the following expression for the electromagnetic 
current: 



cm 



Gn K» Gr 



(29) 



* IT" 



* = 1 



(23) 



The electromagnetic current matrix is the residue of 
the gauged scattering matrix T M at the bound-state pole: 



-M 



(P 2 



* s ; 



M 2 ) ' 



(24) 



where = ^(p i ,q l ,P i ) and = <S>(pf,q f ,Pf) are in- 
and outgoing amplitudes with different momentum de- 
pendencies. On the other hand, Eq. (21) yields 
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1AM 
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M 2 



T 



P 2 + M 2 



(25) 



and by comparing these two equations one obtains the 
electromagnetic current as the gauged inverse scattering 
matrix element between the onshell bound-state ampli- 
tudes: 



(26) 



Electromagnetic current conservation, 



follows from a Ward-Takahashi identity for the five-point 
function (T _1 ) M upon exploiting the onshell condition 
T- 1 * = from Eq. (22). In the limit Q 2 -> 0, Eq. (26) 
reproduces the canonical normalization condition, and 
charge conservation G P E (Q) — 1 is automatically satisfied 
if the Faddeev amplitude is canonically normalized. 

To extract the specific form of the current in the three- 
body approach, one derives from Eq. (20): 



T 



= (K-y -G%= -K- 1 K* K- 1 - G% , 



(28) 



which, upon sandwiching between onshell nucleon ampli- 
tudes and implementing the bound-state relation (22), 



The kernel K = K Gq that appears in the bound- 
state equation (1) of a baryon described by three valence 
quarks is given by 



o=l 



hence the corresponding gauged kernel reads 



K" = Kt 



3 

£ 

a=l 



(k»St\ 
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) 



(30) 



(31) 



where is the dressed quark-photon vertex acting on 

quark (a). In conjunction with Eq. (29) this describes the 
general expression for a baryon's electromagnetic current 
which is displayed in Fig. 6. It comprises an impulse- 
approximation diagram, a component involving the two- 
quark kernel, and two diagrams where the photon couples 
to the qq and the qqq kernels, respectively. In addition 
one must take into account the diagrams with permuted 
quark lines. 

The close relation between the three-body bound-state 
equation (22) and the electromagnetic current (29) makes 
clear that a given truncation must be respected by both. 
Neglecting the three-body irreducible contribution and 
employing a rainbow-ladder two-body kernel which is in- 
(27) dependent of the total momentum, hence K'f = 0, fi- 



nally yields the electromagnetic current 



-ljjr f 



Go 



Go 



V l (32) 



which is depicted in Fig. 7 and employed in the 
present calculation. It is now reduced to the impulse- 
approximation diagram and a further piece that in- 
volves the rainbow-ladder kernel. These contributions 
are worked out in detail in App. D2. The only addi- 
tional component which appears in Fig. 7 is the dressed 
quark-photon vertex. It is computed consistently from its 
rainbow-ladder truncated inhomogeneous Bethe-Salpeter 
equation and described in App. D3. 





+ 



FIG. 7: (Color online) The two types of diagrams which con- 
tribute to the nucleon's three-body current in rainbow-ladder 
truncation, Eq. (32). 



B. Extraction of the form factors 

By virtue of Eq. (18), the electric and magnetic Sachs 
form factors are extracted from the electromagnetic cur- 
rent J M through the Dirac traces 



G E (Q 2 ) 



1 



Tr {J^}, 



(33) 



G M {Q 2 ) = -Tv{J p 1 i r }, 



where r = Q 2 /(4M 2 ), P is the average momentum of 
incoming and outgoing nucleon, and 7^ is transverse with 
respect to P. 

Only the diagrams illustrated in Fig. 7, where the pho- 
ton couples to the respective upper quark lines, need 
to be computed explicitly; the remaining permuted di- 
agrams can be obtained from the relation (D14). The 
incoming and outgoing nucleon amplitudes are given by 
the respective combinations of mixed-antisymmetric and 
mixed-symmetric Dirac amplitudes ^1 and *S>2 that ap- 
pear in Eq. (10). Upon working out the flavor traces one 
finds that only diagrams contribute whose amplitudes in 
the incoming and outgoing state share the same symme- 
try. Thus, both electric and magnetic form factors can 
be decomposed as 



G = G 



AA 



G ss = G 



AA 
IMP 



G 



AA 



^IMP 



(34) 



where G AA denotes a form factor obtained from two 
mixed-antisymmetric amplitudes and G the respective 
contribution with two mixed-symmetric amplitudes. 

The structure of the basis elements in Eqs. (13-14) 
suggests to identify G AA as the predominantly scalar- 
diquark component and G ss as the axial-vector diquark 
contribution to the form factors. While this can serve in- 
tuition, there is however no one-to-one mapping between 
i&x, ^2 and the scalar- and axial- vector diquark contri- 
butions to the nucleon amplitude in the quark-diquark 
model. Ultimately, "Fi and ^2 carry contributions from 
all diquark channels. 

Summing all permuted diagrams, cf. Eq. (D15), yields 
the proton and neutron form factors as the combinations 



2 G AA G^^ G AA 



(35) 



where we again suppressed the subscripts E and M for 
the electric and magnetic contributions. Similarly, one 



obtains the isoscalar and isovector combinations 
G s = G P + G" = G AA + G ss , 



G v = G P - G n = 3G 



AA 



G 



ss 



(36) 



Finally, the flavor-separated up-/down-quark contribu- 
tions in the proton which, owing to charge symmetry, 
equal the down- /up- quark contributions in the neutron, 
are given by 



G u = 2 G p + G n = 3 G AA + G ss . 
G d = G p + 2 G n = 2 G ss . 



(37) 



It is apparent that charge conservation G E (0) = 1, 
G™(0) = requires GjA{0) = Gf s (0) = 1/2, and a 
neutron electric form factor G E (Q 2 ) > can only oc- 
cur if G E A < G E S . Similarly, the neutron's Dirac and 

Pauli form factors obtained from Eq. (19) are negative if 
f aa > F ss 

We further comment upon the relative importance of 
the two diagrams in Fig. 7. From the kinematics in the 
form-factor diagrams, Eqs. (D8-D9), one infers that only 
the relative momenta Pi, Pf of the incoming and outgoing 
amplitudes and the quark momenta pi, p2, p 3 carry a de- 
pendence on the photon momentum transfer Q 2 (in the 
case of pi and P2, through the average momentum P). 
The relative momenta qf as well as the gluon momen- 
tum k in the second diagram are Q 2 — independent, hence 
both diagrams should follow a similar Q 2 — evolution. In- 
deed we find that in the accessible photon-momentum 
range Q 2 < 7 GeV 2 the ratios 



R{Q 2 ) := -G k (Q 2 )/Gtmp(Q 2 



(38) 



tors, R p A f 
factor, R P E 



depend only weakly on Q 2 : for the magnetic form fac- 
0.3 . . . 0.5 and for the electric proton form 
' 0.4...1. For the electric neutron form 
factor charge conservation implies Gk(0) = — Gtmp(0), 
hence R E (0) = 1. While the kernel contribution to G E 
stays positive at all Q 2 > 0, the impulse-approximation 
diagram is negative at Q 2 = and becomes positive at 
Q 2 « 1 GeV 2 above which it dominates the form fac- 
tor. Both contributions in Fig. 7 are therefore hard, and 
the diagram which involves the qq kernel and necessarily 
augments the impulse-approximation diagram to ensure 
current conservation stays important at large momentum 
transfer. 



IV. RESULTS AND DISCUSSION 

We computed the electromagnetic form factors of the 
nucleon within the decomposition of Fig. 7 for the electro- 
magnetic current. We emphasize that the present setup 
allows to determine all ingredients of the electromag- 
netic current without any further approximations: the 
rainbow-ladder kernel that appears in its definition is 
given in Eq. (8); the dressed quark propagator is obtained 
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Q 2 /M 2 



Q 2 /M 2 



FIG. 8: (Color online) Electromagnetic form factors of the nucleon as a function of the photon momentum transfer Q 2 . 
The bands indicate a variation of r\ £ [1.8,2.0]. The experimental data for and Glf are from Refs. [83-91] and [92-97], 
respectively. The references for G P E are given in Fig. 10 and those for G V M can be found in [1]. 





Mat 






K V 




' E 




r p 

M 


r n M 


Faddeev 
Exp. 


0.94 
0.94 


2.21(1) 
2.79 


-1.33(1) 
-1.91 


2.54(2) 
3.70 


-0.12(1) 
-0.12 


0.75(3) 
0.89 


-0.01 
-0.12 


0.72(2) 
0.86 


0.72(2) 
0.87 



TABLE II: Results for the nucleon mass and static electromagnetic properties compared to experiment. The parentheses 
indicate the dependence on the infrared parameter r\ of Eq. (9). Mn is given in GeV, the magnetic moments are expressed in 
nuclear magnetons, and the charge radii are given in fm. 



from its Dyson-Schwinger equation (7); the nucleon am- 
plitude is the solution of the covariant Faddeev equa- 
tion (2); and the quark-photon vertex is obtained consis- 
tently from its inhomogeneous Bethe-Salpeter equation, 
also solved within a rainbow-ladder truncation, and is de- 
scribed in App. D 3. Apart from the current-quark mass, 
the single parameter that enters the equations and con- 
trols the results is the infrared scale A which appears in 
Eq. (9) in connection with the rainbow-ladder trunca- 
tion, where it was fixed to reproduce the experimental 
pion decay constant. In the chiral limit, A is the only 



relevant scale in the system, and all mass-dimensionful 
quantities scale with A. This entails that dimensionless 
form factors become independent of this scale. One can 
further investigate the impact of the infrared properties, 
described by the width parameter rj (cf. Fig. 3), on re- 
sulting observables. This is indicated by the bands in 
Figs. (8-10), and in the same way as n— and p— ground 
state properties are not very sensitive to r), this sensitiv- 
ity is found to be weak in the case of nucleon form factors 
as well. 
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Low-momentum behavior 



Eq. (39) with the following reference mass [50] : 



Fig. 8 shows the results for the nucleon electromagnetic 
form factors, calculated at the physical u/d point corre- 
sponding to m v = 140 MeV and compared to experiment. 
We find a remarkable agreement with the experimental 
data above a photon momentum transfer Q 2 > 2 GcV 2 . 
This is the region where pseudoscalar-meson cloud ef- 
fects should vanish as the photon probes the nucleon at 
length scales much smaller than the typical size of pionic 
correlations and thereby essentially reveals the nucleon's 
quark core. An illustrative example is the neutron's elec- 
tric form factor where our curve is close to the new Hall- 
A data from JLAB [91], whereas the bump at low Q 2 
is completely absent in our result. This suggests that 
the low-O/ 2 structure of G% predominantly owes to vir- 
tual pion-cloud components. Similar effects are found in 
the proton's and neutron's magnetic form factors: below 
Q 2 ~ 2 GeV 2 , the results underestimate the data and, at 
vanishing photon momentum, yield magnetic moments 
that are 20 — 30% smaller than the experimental val- 
ues, see Table II: for the proton and neutron, they read 
H p = 2.21(1) and // l = -1.33(1), where the brackets 
denote the sensitivity to r\. 

Fig. 9 displays the current-quark mass evolution of our 
calculated nucleon magnetic moments and charge radii, 
again expressed in terms of the squared pion mass which 
is obtained from its Bethe-Salpeter equation. The upper 
left panel depicts the dimensionless magnetic moments 
fjP' n = G^ n (0) which exhibit a moderately rising be- 
havior with the pion mass. In the upper right panel of 
Fig. 9 the quark-mass dependence of the isovector and 
isoscalar anomalous magnetic moments k v = n p — K n 
and k s = k p + K n , expressed in static nuclear mag- 
netons, is compared to recent lattice data. A very 
similar current-mass dependence emerges in the quark- 
diquark model [50]; however, compared to those results, 
the gluon-interaction kernel in the explicit three-quark 
framework leads to an overall reduction of 10 — 15% for 
the magnetic moments. 

For a comparison of magnetic moments in units of 
static nuclear magnetons with experiment or lattice, one 
must bear in mind that our calculated nucleon mass dif- 
fers (although marginally, cf. Fig. 5) from that obtained 
on the lattice at higher pion masses. The nucleon's mag- 
netic moments are given by 



Mdim 



2M 



G P M n (0) 



2AL- 



r<P' n (r\\ ^ cx p 



(39) 



where M is running with the current-quark mass; hence 
their values in static nuclear magnetons are given by 
the bracket in Eq. (39). The unambiguous comparison 
of magnetic moments is that of the dimensionless value 
Gm (0) , whereas the corresponding magnetic moments in 
static nuclear magnetons will compare differently if the 
nucleon masses in both approaches do not coincide. To 
account for this, we plot k v,s in Fig. 9 by replacing M in 



2 cf (m 2 ) = Af 2 +^) 2 (l + /(m 2 )), (40) 



M 



where /(to 2 ) = 0.77/(1 + (to^/0.65 GeV) 4 ) and M = 
0.9 GeV. Eq. (40) reproduces the experimental nucleon 
mass at m, = 0.14 GeV, approaches the heavy-quark 
limit via M — > 3to„-/2 and describes the dynamical lattice 
results for the nucleon mass in Fig. 5 reasonably well. 

Chiral cloud corrections to core magnetic moments can 
be estimated from the pion-loop contributions 5fi p ' n = 
lJ>,n _ ^p^ c in heavy-baryon chiral effective field theory. 
In that framework the sum of the two diagrams where 
the photon strikes a pion in the nucleon and leaves an 
intermediate nucleon or A— baryon as a spectator is given 
by [101, 102] 



Sfi p ' n =±X dx — u{x) 



C 2 w(2w + l) 
9ffl (^ + 1) 2 



, (41) 



where x = |fe|/(AM) and u> — \J x 2 + to 2 / (AM) 2 are 
the pion momentum and intermediate pion energy, nor- 
malized by the physical N-A mass splitting AM = 
0.29 GeV. The constants are given by C = -2 -(0.76) 
and A = g 2 A M N (AM)/{3Tr 2 f 2 ) with /„ = 131 MeV 
and gA = 1.26. Choosing an identical dipole ansatz 
u(x) = [1 + (|fc|/A) 2 ]~ 2 for the NNn and VAtt vertex 
dressings with a regulator A = 0.8 GeV yields the pion- 
loop contribution Sfi p,n — ±0.61 at the physical pion 
mass; however, the inclusion of further meson-loop dia- 
grams can diminish this value [102]. Upon subtracting 
this result from the experimental values one obtains the 
core estimates /j% oic = 2.18 and /x™ orc = —1.30. This is 
remarkably close to our calculated results /i p = 2.21(1) 
and n n = -1.33(1). 

Since the simplest pion-loop contributions to proton 
and neutron magnetic moments of Eq. (41) carry an op- 
posite sign, their total cancels in the isoscalar combina- 
tion. This implies k s w K tnre-> i- e - isoscalar magnetic 
moments should be roughly undisturbed by pion-cloud 
corrections. The experimental isoscalar anomalous mag- 
netic moment is small and negative, k s — —0.12, and re- 
produced by the Faddeev calculation, cf. Table II, with 
only a small model dependence. 

The electromagnetic radius corresponding to a form 
factor G{(Q 2 ) is defined as 



6 dG t 
G~(0j dQ 2 



(he) 



(42) 



Q 2 =0 



where he = 0.197 GeV fm, and Gi(0) in the denomina- 
tor is dropped for the neutron's form factors and F™ 
as they vanish at Q 2 = 0. The lower panels in Fig. 9 
show the squared isovector Dirac and Pauli radii (r^) 2 
and (r^) 2 , compared to lattice data. Once again, the 
unambiguous comparison with the lattice is that of the 
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FIG. 9: (Color online) Static electromagnetic properties of the nucleon as a function of the squared pion mass, for 77 = 1.8±0.2. 
Upper left panel: dimensionless magnetic moments of proton and neutron. Upper right panel: isovector and isoscalar magnetic 
moments in units of static nuclear magnetons. Lower panels: squared isovector Dirac and Pauli radii, given in fm 2 . The lattice 
data are from Refs. [72-76, 98]. Stars denote experimental values [99, 100] at the u/d— quark mass whose position is indicated 
by dashed vertical lines. 



dimensionless values r?M 2 , hence we rescale our calcu- 
lated radii via 

2 2 ( M calc \ 

' resc ' calc I n r J 5 v^ / 

with MR e f from Eq. (40). The plot shows a satisfactory 
agreement with the lattice data at larger pion masses 
where the pion-cloud dressing effects of the nucleon are 
diminished. Compared to the reduction of the magnetic 
moments, the isovector radii from the Faddeev approach 
are virtually identical to those obtained in the quark- 
diquark model [50] . The absence of a meson cloud is fur- 
ther signaled by missing chiral curvature as the charge 
radii of the nucleon, surrounded by a pion cloud, would 
diverge in the chiral limit. At the u/d mass, the ex- 
perimental radii rL, r p M and r 1 ^ are underestimated by 
15 - 20%, cf. Table II. 



To summarize this section: the combined behavior of 
the nucleon's magnetic moments and the electromagnetic 
radii, compared to experiment and lattice data, provides 
strong evidence that the main missing ingredients to nu- 
cleon form factors in a rainbow-ladder calculation are 
those induced by chiral-cloud corrections. Thus, Fig. 9 
outlines the characteristic features of a current-mass de- 
pendent nucleon quark core. 



B. Large-momentum behavior 

With the body of new experimental data at higher pho- 
ton momentum transfer, the large-Q 2 behavior of the nu- 
cleon's electromagnetic form factors has attracted consid- 
erable attention in the past decade. Perturbative QCD 
predicts the scaling behavior of the Dirac and Pauli form 
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Q 2 /M 2 Q 2 /M 2 

FIG. 10: (Color online) Q 2 — evolution of nucleon electromagnetic form factor ratios. The bands denote the n variation of 
Fig. 8. Upper left panel: Electric proton form factor normalized by the dipole (see text) and compared to experimental data 
from Refs. [103-109]. Upper right panel: Pauli-to-Dirac ratio for the up-quark, normalized by the experimental value k u = 1.67, 
with data from Ref. [110]. Lower panels: weighted Pauli-to-Dirac ratios of Eq. (46) for proton and neutron. The experimental 
data for are the same as in the upper left panel and those for i?2i are identical to G\ in Fig. 8. 



factors as [111] 

F 1 ^l/Q\ F 2 ^l/Q 6 , Q 2 F 2 /F 1 ^ const. (44) 

Correspondingly, the Sachs form factors would scale as 
Ge,Gm — > 1/Q 4 which implies that the ratio Ge/Gm 
becomes constant. These predictions have come under 
scrutiny with the polarization-transfer measurements, 
where the ratio G E /G P M shows roughly a linear decrease 
with Q 2 and points toward a zero crossing at some larger 
value of Q 2 . 

The discrepancy between the perturbative prediction 
and the experimental data has been attributed to he- 
licity non-conservation and the presence of non-zero 
quark orbital angular-momentum in the nucleon ampli- 
tude [10, 112, 113]. Indeed, with an updated perturbative 



prediction for F 2 that accounts for wave-function compo- 
nents with orbital angular momentum [9], 

Q 2 {F 2 /F x )/\n 2 (Q 2 /A 2 ) -> const., (45) 

the onset of perturbative scaling in the proton's Pauli to 
Dirac ratio appears to happen already at comparatively 
low photon momenta, cf. Fig. 10. 

While the large-Q 2 behavior is in principle accessible 
in the Faddeev approach, its current implementation is 
limited by an upper value of the photon momentum. The 
quark propagator obtained from the DSE inevitably de- 
velops a certain singularity structure which can be char- 
acterized by a pole mass m q [27] . The photon probes the 
quark propagator S(p) within a parabola in the complex 
p 2 — domain that grows with Q 2 and is thereby limited by 
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the value of m q . In absence of a sophisticated method 
to account for these singularities, a bound-state ampli- 
tude can be computed from the Faddeev equation only if 
M < 3m q , and the form- factor calculation is restricted to 
Q 2 < 4 ((3m 9 ) 2 - M 2 ) fa 7 ... 9 GeV 2 , depending on the 
infrared parameter r\. In proximity to this upper bound 
the results become sensitive to the numerics and require 
an increasingly better resolution of the Faddeev ampli- 
tude. Thus, with presently available accuracy, the results 
in this region should be interpreted with caution. 

The upper left panel of Fig. 10 shows the electric form 
factor of the proton normalized by the dipole ansatz 
G D = (l+c D Q 2 /M 2 )- 2 , where c D = M 2 xp /(0.71 GeV 2 ). 
It exhibits a fall-off in Q 2 which signals the presence of 
orbital angular momentum in the Faddeev amplitude, 
cf. Table I. While we do not yet deem our numerical re- 
sults fully reliable in the large-momentum region, Fig. 10 
demonstrates that a zero crossing can occur quite natu- 
rally within the Faddeev approach. In fact, if Eq. (45) 
is valid at already moderate photon momenta it eas- 
ily accommodates a zero crossing for G P E /G P M ; namely, 
through Eq. (19), at ln 2 (Q 2 /A 2 ) = AM 2 /const. Quark- 
diquark model studies typically find a zero crossing as 
well, where its location depends on the model parame- 
ters in the calculation [47, 113]. 

In connection with this, an essential component to en- 
sure the Q 2_ev °hition in Figs. (8) and (10) is the quark- 
photon vertex solution from its inhomogeneous BSE. 
While the Ball-Chiu ansatz of Eq. (D19) alone is suffi- 
cient to guarantee current conservation, it misses impor- 
tant transverse parts in the vertex which are dominated 
by a vector-meson pole in the vertex at Q 2 = —m 2 [114]. 
It has been noted in several works that the Ball-Chiu ver- 
tex contributes only ~ 50% to squared hadronic charge 
radii [50, 114-116]. That statement stays true here. 
While the magnetic moments are not very sensitive to 
these transverse contributions, the Q 2 — evolution changes 
dramatically: at Q 2 /M 2 — 3, the form factors G E , G P M 
and G^j obtained from the Ball-Chiu ansatz alone over- 
estimate the results in Fig. 8 by a factor of 3, whereas 
G E develops a zero crossing at Q 2 /M 2 ps 4. This makes 
clear that the position of a zero crossing might well de- 
pend on the truncation, as effects beyond rainbow-ladder 
could certainly impact upon the transverse structures of 
the quark-photon vertex and thereby play a delicate role 
in the form factors at larger Q 2 . 

The lower panels of Fig. 10 display the weighted ratio 
of Pauli to Dirac form factors for proton and neutron, 

_ Q 2 F 2 1 , . 

21 ' M 2 F^^icRQ 2 ^! 2 )' K ' 

where c R = (M oxp /(0.3 GeV)) 2 . Eq. (45) entails that 
the ratios -Rfi™ should approach constant values at large 
photon momenta. While such a near constancy in R^i 
is observed already at moderate Q 2 , the neutron's ratio 
does not yet follow this trend. Irrespective of that, we 
find a fair agreement between our results and the data. 



An intriguing feature which was recently noted [110] 
is the behavior of the up- and down-quark ratios F% j 
and F^/Fi above Q 2 ~ 1 GeV 2 : compared to their pro- 
ton and neutron counterparts, they are roughly constant, 
whereas F% / Ff rises sharply below that value toward its 
static value i/a, as shown in the upper right panel of 
Fig. 10. Combined with the suppression of the d— quark 
form factors compared to those of the u— quark this ex- 
plains the observed behavior of an d -R21 ■ The decline 
of Ff j F" and F$ / F% is responsible for the negative form 
factors F™ 2 of the neutron. In the infinite-momentum 
frame it is related to a concentration of the up-quark dis- 
tribution in the proton, and that of the d— quark in the 
neutron, in impact-parameter space [15]. In the quark- 
diquark model this suppression owes to the singly rep- 
resented d— quark in the proton that mainly interacts 
with the photon in combination with an axial-vector di- 
quark in a uu configuration [47]. The upper right panel 
of Fig. 10 shows that the constancy of F% /F™ also main- 
tains for higher values of the photon momentum, whereas 
the rise at small Q 2 is practically absent and thereby 
likely to be related to pion-cloud effects. 

We finally remark that there is no qualitative differ- 
ence in the form factors' Q 2 — evolution up to the strange- 
quark mass: both the form factor curves of Fig. 8 as well 
as their ratios in Fig. 10 essentially retain their shape 
when plotted as a function of the dimensionless variable 
Q 2 /M 2 . This feature was also observed in Ref. [51] for 
the A electromagnetic form factors in the quark-diquark 
context. It can be understood from the correlations in the 
Faddeev amplitude, cf. Table I: between the up-/down 
and strange-quark regime, orbital angular momentum 
contributes one-third to the nucleon spin and this con- 
tribution only slowly decreases with rising current-quark 
mass. For the quark masses investigated here, the dis- 
appearance of p— wave components in the nucleon's rest 
frame which is expected to have a sizeable impact on 
nucleon form factors is therefore not yet visible in our 
results. 



V. CONCLUSIONS 

We presented a calculation of nucleon electromagnetic 
form factors in the Poincare-covariant Faddeev approach, 
based on the Dyson-Schwinger equations of QCD. We 
employed a rainbow-ladder truncation which explains the 
binding in the nucleon through iterated dressed gluon ex- 
change between the quarks. Thereby all ingredients of 
the equations are completely fixed and no further model 
assumptions need to be made. The single parameter of 
the approach is a scale which is fixed to reproduce the 
pion decay constant. The nucleon's electromagnetic cur- 
rent in terms of three interacting quarks was derived. 
Since pion-cloud effects are not implemented, our results 
represent the nucleon's quark core. 

We find a good agreement between our form-factor re- 
sults and the experimental data above Q 2 pa 2 GeV 2 . 



14 



The form- factor ratio G P E /G P M falls off with larger Q 2 and 
the results for the Pauli to Dirac ratios closely follow the 
data. These features can be attributed to a significant 
amount of orbital angular momentum which appears in 
the solution for the nucleon's amplitude and persists up 
to the strange quark mass, and most likely much farther 
beyond. In the low-Q 2 region clear signals of missing 
pion-cloud effects are present: the charge radii and mag- 
netic moments underestimate the data, and the enhanced 
low-Q 2 structure in the neutron's electric form factor is 
absent. At larger quark masses, where pionic effects do 
no longer contribute, our results for the nucleon mag- 
netic moments and electric charge radii are comparable 
to those obtained in lattice QCD. 

The present framework does not include genuine three- 
quark interactions. However, our results suggest that 
such contributions are indeed small and that it is es- 
sentially the quark-quark correlations which contribute 
most of the nucleon's binding. The overall agreement 
between our results and those obtained in the quark- 
diquark model provides further evidence for the quark- 
diquark structure of the nucleon, and it implies that 
scalar and axial-vector diquark degrees of freedom can 
account for most of its characteristic features. 

In view of a better understanding of the chiral and low- 
Q 2 structure of baryons, our approach must be improved 
by incorporating missing chiral cloud effects. This is a 
challenging task as it necessitates a consistent descrip- 
tion of baryons beyond the rainbow-ladder truncation. 
Steps in that direction are planned. Moreover, an inves- 
tigation of nucleon electromagnetic form factors at large 
photon momentum as well as nucleon-to-resonance tran- 
sition form factors is desirable. The results presented 
herein are encouraging, and we are confident that the co- 
variant Faddeev framework provides a suitable starting 
point for future studies of baryon structure. 
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Appendix A: Conventions and formulas 

1. Euclidean conventions 

We work in Euclidean momentum space with the fol- 
lowing conventions: 

4 

P-1 = ^2PkQk, p 2 =p-p, i> = vi- (Al) 
k=i 

A vector p is spacelike if p 2 > and timelike if p 2 < 0. 
The hermitian 7— matrices 7^ = (7^ satisfy the anti- 
commutation relations {7^,7"} = 2 5^, and we define 

v* v = ~kf,'f], 7 5 = -7V7V • (A2) 
In the standard representation one has: 

k= ( -io*\ 4= (t \ Wo 1\ 

where a k are the three Pauli matrices. The charge con- 
jugation matrix is given by 

C = 7V, C T = C f = C- 1 = -C , (A3) 

and the charge conjugates for (pseudo-) scalar, (axial-) 
vector and tensor amplitudes are defined as 

r(p,P) :=CT(-p,-P) T C T , 
T a (p,P) :=-CT a (-p,-P) T C T , (A4) 
T 0a (p,P) :=CT afi (-p,-P) T C T , 

where T denotes a Dirac transpose. Four-momenta are 
conveniently expressed through hyperspherical coordi- 
nates: 




V z J 
and a four-momentum integration reads: 

OO 1 1 27T 

/ - (2 ^2/* 2p2 /^ v/l ^ 2 /^/ # - 

p 0-1 -10 

2. Transverse identities 

The Faddeev- equation formalism simplifies when work- 
ing with quantities that are transverse to the nucleon 
momentum P, where P 2 = —M 2 . Covariant relations 
between transverse objects bear a close resemblance to 
those in three-dimensional Euclidean space: in the nu- 
cleon's rest frame one has P — e^, and quantities trans- 
verse to e4 reduce to the usual three-dimensional form. 
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In a form factor calculation, typically carried out in the 
Breit frame, the average momentum (Pi + P/)/2 of in- 
coming and outgoing nucleon is at rest. We define the 
transverse projector with respect to P: 



= — p^p v 



(A6) 



as well as 7— matrices transverse to P: 

7 £ = T£V = 7 M - P"? , (A7) 

which satisfy {7^,7^} = 2Tp V and 7^7^ = 3. Then 7^, 
75 and f pairwise anticommute. Similarly, we denote 
momenta transverse to P by 



p»=T£ v p v = p^-p-PP» 



(A8) 



and define 



[ 7 £, 7 £]=i(Tjf - 7 £ 7 £). (A9) 



We make frequent use of the positive- and negative en- 
ergy projectors 

A W (P) = (l+wf)/2, w = ± (A10) 
which obey the relations 

A u 75 = 75 A_ w , 



K 1t = 7 £A_ 



(All) 



It is convenient to define a transverse e— tensor by 

a&l c aB^\ pA 

t T .— t r , 

where e 1234 = 1 which satisfies 



(A12) 



flap fj,pcr map rpficr rpacr rpfip 

Crp €.rp - — _L p P P P ' 

(1VOL cy rpQ.(3 

£rp €-rp Zi _£ p , 



(A13) 



and 



75 £t j i-- l p1t + 1 p1t- 1 p1t 



where the antisymmetric color part is normalized to 1. 
(^n) a pys(p,q, P) and (F n ) abcrf are the spin-momentum 
and flavor amplitudes which transform as doublets under 
the permutation group S 3 , with entries for n = 1,2. The 
Dirac amplitudes ^ n carry 3 spinor indices for the quark 
legs (a, j3, 7) and one spinor index S for the nucleon. They 
are mixed-antisymmetric (^i) or mixed-symmetric (^2) 
under exchange of the indices a, /3 and corresponding 
quark momenta. Likewise, the two isospin- 1 ^ flavor ten- 
sors F„ carry 3 isospin indices (a, b, c) for the quarks and 
one (d) for the nucleon and are mixed-(anti-)symmetric 
in a, b. They read: 



1 . 



F 2 = -V6 <TJff2 



(B2) 



where Oi are the Pauli matrices and the tensor product 
is understood as 



(f ' ®g)abcd 
(A ®/2)(5l® 52 ) 



fab Qcd 1 

(/lffl) ® (/2 52) 



(B3) 



They are normalized to (P n ,) bad* c(F n ) abed = <W<W- To 
project onto proton or neutron flavor states, the index d 
is contracted with either of the two isospin vectors (1,0) 
or (0,1). 



1. Permutation-group properties 

The Pauli principle requires the full Dirac-flavor-color 
amplitude of Eq. (Bl) to be antisymmetric under ex- 
change of any two quark legs, i.e. under a combined per- 
mutation of Dirac, flavor and color indices and the corre- 
sponding exchange of momenta. The color amplitude in 
Eq. (Bl) is an antisymmetric singlet. The transformation 
properties of the doublets (B2) follow from Fierz identi- 
ties. With the constraint of a fully symmetric Dirac- 
flavor part one thus obtains the transformation behavior 
of the Dirac amplitudes. We abbreviate the 6 possible 
permutations by 



_ a 3 7 
IT 7t IT ' 


(A14) 


(123) : 


F„ = 


i^n^abcd 7 




{^ n ) a 0js(p,q,P) 




(A15) 


(231) : 


K = 


(F njbcad 7 


= 

n 




7£7t75 4"T = 2 7?, 


(A16) 


(312) : 


K = 


{^n)cabd 7 


tfr" - 

n 


{^n) la 8 S (p",q",P) 


It It It 75 f = 6 • 


(A17) 


(213) : 


F„ = 


i^njbacd 3 


$ — 


(Vn)pa-ys(P,-q,P) 






(321) : 


K = 


(Fn)cbarf 7 


= 

^ n 


(^„)^ aS (p',~q',P) 


Appendix B: Nucleon amplitude 




(132) : 


K = 




\j>" - 

^ n 


(V n ) a ^ S (p",-q",P) 



The nucleon's covariant three-quark amplitude includ- 
ing its full Dirac, flavor and color dependence is given 
by 



where for instance (231) denotes a permutation 
{pi,P2,P3} — ► {P2, J>3, Pi}- The resulting permuted rela- 
tive momenta read: 



*(^,P)=[g*nFj^, 



(Bl) 



P 
q' 



-q- f > 

q , 3p 
2^4' 



P =9- 2 ' 

// _ q 3p 



(B4) 



16 



The Dirac and flavor doublets (for the latter: replace \? 
by F) transform under these permutations as 



16 tensors T,j with i,j = 1 . . .4 which are detailed below. 
The Dirac amplitudes are then reconstructed via 



* = M'V = M"^" = = M 1 ^ 1 = M"^" , (B5) 

where the transformation matrices acting upon the index 
n axe given by A4 — diag(— 1, 1) and 



m' = -( ~i 

2 \ V3 -1 



1 / 1 



-V3 
2 I -73 -1 




Because of their orthogonality a scalar product of two 
doublets, such as the bracket in Eq. (Bl), transforms as 
a symmetric singlet under a permutation. 

It is convenient to define Lorentz-invariant momentum 
variables with definite transformation properties under 
the permutation group, cf. Ref. [117], a feature which is 
not manifest in the set {p 2 , q 2 , z , z\, z 2 }. Indeed one can 
construct a symmetric singlet 



? 2 

V Q ^ 



(B6) 



and two further doublets 




^-1 




(B7) 



on the domain (— 1, 1) e R, with £ := (4<j 2 )/(3p 2 ) and 



P • g = 2i2 2 + z Q J 1 - zf W 1 - z. 



(B8) 



Under a permutation of momenta (B4) the doublets again 
transform via Eq. (B5). 



2. Orthonormal basis 

The mixed- (anti-) symmetric Dirac amplitudes *5f n in 
Eq. (Bl) have the general structure of four-fermion Green 
functions with positive parity and positive energy. As 
such they are described by 64 covariant basis elements 
which we denote by 



Xi, 



Tij (A w75 C<g> A+) . 



^2,iju ■= (75 ® 75) Xl^jw 



(B9) 



where A w (w = ±) are the positive- and negative-energy 
projectors from Eq. (A10), and restriction to positive en- 
ergies for the nucleon eliminates the occurrence of A_ on 
the right-hand side of the tensor product. The depen- 
dence on the relative momenta p and q is carried by the 



#„(p,?,P) = fnMjUt)*k,ijUr,s,P) , (BIO) 

kijuj 

where the coefficients f n ,kiju>{t) depend on the 5 Lorentz- 
invariant momentum variables t := {p 2 , q 2 , zq, z±, z 2 }. It 
is convenient to express the momenta {p, q, P} through 
orthogonal unit vectors {r, s, P} which satisfy r 2 = s 2 = 
P 2 = 1 and r ■ s — r ■ P = s ■ P = 0. This is achieved in 
a covariant way via 



■ p T 



P - ZiP 



q T 



q-z 2 P 



which, in the nucleon's rest frame (P = e^j and the co- 
ordinate representation (A5), (for p^\ replace z — > z\ 
and y — > 1; for q^: p 2 — > q 2 , z — > z 2 , y — > z , <fi — > 0) 
assume the simple momentum-independent form r = 
and s = e 2 . 

Let us now define the basic spin structures 



r J= i... 4 = {1, 75^, f, f\ 



(Bll) 



with v 11 = et^ al3 r a s 13 and by virtue of Eq. (A15): 75 j) = 
f $ f. The assignment = T t <E> Tj would yield a com- 
plete but non-orthogonal basis (B9). An orthonormal ba- 
sis which also corresponds to a partial-wave decomposi- 
tion in the nucleon's rest frame was presented in Ref. [34] . 
It reads 



Tij = i <g> r 



J > 



T 2j = -^ 7 ?®7?r j , (B12) 



with Tj from Eq. (Bll), and 



T 



1 



3J ~ V6 



,ry : 



■It' 

where T a , Tf (j = l. 



T 



1 

71 



7?®7^rf : 



(B13) 



. 4) are defined by 



3r°y _ 7 q 

(3 ^-7^)75 
3 r a - 7 £ / 
3 s Q - 7^ 

respectively, with 



and 



cf = (2s a s fi -T^)e 



r a s P + s a r f3 

ad 
Ci 75 

aB 
C 2 75 



1 p 



a/3 



The and L"^ satisfy the relation 



7 £i7 = 7 £ 7 £rf = 0. 



(B14) 



(B15) 
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With this choice, the basis elements in Eq. (B9) obey 
the orthogonality relation 



} 



Skk'Sii'Snj'Su, 



(B16) 



Moreover, they are eigenstates of the total quark spin 
S 2 and orbital angular momentum L 2 with eigenvalues 
s(s + 1) and 1(1 + 1). The s = 1 /2 subspace is spanned 
by the elements in Eq. (B12); the s — 3 / 2 subspace by 
those in Eq. (B13). Tn and T21 are s— waves (/ = 0), 
T31 and Tij with j = 1 ... 4 are d— waves (I = 2) and 
the remaining ones are p— waves (I = 1). We will carry 
out the partial-wave decomposition in detail in the next 
subsection. 

We note that the Tjj of Eqs. (B12-B13) also consti- 
tute a complete orthonormal basis of a four-quark Green 
function with positive parity, however with Eq. (B9) gen- 
eralized to the form 



X-2,ijuu' ■= (75 ® 75) Xl,^, 



(B17) 



where we dropped 75 C to indicate the case of a quark- 
antiquark four-point function. 



3. Partial-wave decomposition 

The Dirac basis elements X^y^ of Eq. (B9) can be clas- 
sified with respect to their quark-spin and orbital angu- 
lar momentum content in the nucleon's rest frame which 
we will demonstrate explicitly in this appendix. Only 
the total angular momentum j = 1/2 of the nucleon is 
Poincare-invariant while the interpretation in terms of 
total quark spin and orbital angular momentum will dif- 
fer in every frame. The spin is described by the Pauli- 
Lubanski operator: 



1 



eti a0 J o 



(B18) 



where we chose the total momentum P to be normalized, 
cf. Eq. (A12). As generators of the Poincare algebra, J^ v 
and P M satisfy the usual commutation relations. The 
eigenvalues of the square of the Pauli-Lubanski operator, 



W 2 



1 



rpfia rpV/3 JftV J 



a/3 



iO' + i) 



(B19) 



define the spin j of the particle. For a system of three 
particles with total momentum P and relative momenta 
p and q, the total angular momentum operator consists 
of the total quark spin S and the relative orbital angular 

Upon subsuming them into 



momentum L 



'(9) 



Lorentz-covariant operators 



5" 
re 



: • , (a af} 0101 + perm. 
§4^(p a d£-/6£)l<g>lC 



(B20) 



_ i ua/3 
2 fc T 



(q a d$ - q d*) 10101 



L^ q y one can show that the 



with 

orthonormal basis elements of Eq. (B9) are eigenstates 
of the operators 

S 2 



1 ® 1 <g> 1 + \ (a% v ® a^ v ® 1 + perm.) , (B21) 



and L 2 = L 2 p) + L 2 q) + 2 L (p) -L( q ), where 



T 2 



2p T -d p + (p£ p v T - p% T£") dt d; 



L 2 q) =2q T -d q + («# <# - q 2 Tjf) d% 8" , (B22) 

L (p)' L (q) = (Pt Qt~Pt- qr Tp V ) dp d q , 

with eigenvalues s(s + 1) and 1(1 + 1). 

In the nucleon's rest frame, where P = the opera- 
tors (B21-B22) assume the meaning of total quark spin 
and orbital angular momentum. Here S 2 reads explicitly 



1 + perm.) (B23) 



where Ej = -^£ijk&jk- Choosing hyperspherical coordi- 
nates {p 2 , z, y, (f)} for p^ and {q 2 , 2;', y', 0'} for g M ac- 
cording to Eq. (A5), defining p, q as the three- vectors 
corresponding to px and qx, e.g. 



y z sm <; 

P = Vp 2 Vl-z 2 ( v/1 - y 2 cos< 

y 



(B24) 



and expressing the gradient and Laplacian in spherical 
coordinates, the operator L 2 is the combination of 

L 2 P ) = 2p • V p + p h (p ■ V p ) V* - p 2 A p 



2 2/ a 2/ -(i- 2/ 2 )a 2 



3 
1-y 2 



L 2 „^ = 2q ■ V q + q k (q ■ V q ) V* - q 2 A,; 



= 2y '^-(i-^)s 2 ,- i ^ 



and 



l (p) -hi) = p k (i ■ v p) V* - (p ■ g)(V p • V,) = 

= - cos(^ — 4>')\f I — y 2 \J 1 — y' 2 d y 8y> 
yy' cos(0 - 0') 



sinfii — 6') 



y'(l-y 2 )d v dt,-y(l-y' 2 )d y ,dt 



While these expressions are expedient, e.g., for use 
with a symbolic programming language, Eqs. (B21-B22) 
allow for a partial-wave decomposition in a more direct 
way. To find the eigenstates of S 2 one substitutes in 
Eq. (B21) the relation 



-2 7^75 



7t 7e 



f (B25) 
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that follows from Eq. (A15). Then, with the anticommu- 

tation properties of 7^, 75 and f and repeated use of 
Eq. (A15) it is straightforward to show that 

with s = Y 2 for i = 1, 2 and s = 3 /2 for i = 3, 4. The proof 
only relies on the generic form (B12-B13) of the T^, 
together with the property (B15), but not the specific 
form of the matrices Tj, and T"^. 

The eigenstates of L 2 can be determined independently 
of the Dirac structure since the operators L 2 p ^ , L 2 q ~. and 
L( p ) ■ only act upon the relative-momentum depen- 
dence. The elements Tn and T21 are independent of the 
relative momenta, hence they carry orbital angular mo- 
mentum / = 0. All other basis elements in (B12-B13) 
can be expressed through contractions of 



r a , s a , r a s^ s a r^ r a , s a , 
r^s 7 , s a s fi r 1 



(B26) 



L 2 r a 


= 2r a , 




L 2 s a 


= 2s a , 




L 2 r a r< 3 


= 6rV 


-2Tf 


L 2 s a S P 


= 6s a s P 


-ITf 


L 2 r a s? 


= 4r a / 


+ 2sV 



with appropriate Dirac structures, for instance: T13 = 
(1 ® 7^) r a . Applying the operator L 2 on (B26) yields 
for example 



(B27) 



and the results for higher powers of momenta follow from 
the relation 

L 2 (fgh) = f L 2 (gh) + g L 2 (fh) + h L 2 (fg) 
-(gh)L 2 f-(fh)L 2 g~(fg)L 2 h. 
From (B27) one constructs the eigenfunctions for I = 1: 

(B29) 

The unit vectors r a , s a and v a appear in the basis ele- 
ments Tij for i = 1,2,3 and j = 2,3,4 which therefore 
carry 1 = 1. For I = 2 one obtains the eigenfunctions 

r^ + sV, 3r a r^-Tf, 3s a s p -Tf (B30) 

that appear in T31 , T41 and T42 , and also the eigenfunc- 
tions 

A m (2 8 a sP - if) r 7 = -<S? etf 7 , (B31) 



A lM {2r^-Tf)s^ = -cr 



2 C T 



(B32) 

which occur in T43 and T44. Here A[^] denotes an anti- 



a0 



symmetrizer with respect to the indices /3, 7 and the c 12 
are defined in Eq. (B14). A second set of eigenfunctions 
with A\ a -A instead of A\p-A yields no new information. 

To summarize, the 16 elements Ty can be categorized 
into two s— waves, nine p— waves and five d— waves. In 
the final basis elements of Eq. (B9) these numbers are 
multiplied by four through the additional indices k = 1, 2 
and lu = ±. 



Appendix C: Faddeev equation: updated solution 
strategy 

In this appendix we detail an updated solution strat- 
egy for the Faddeev equation which considerably reduces 
the numerical effort. We start by writing the Faddeev 
equation for the Dirac part of the nucleon amplitude as 



$C0 + *(2) + $(3) 

n. 1 I n 1 n. 



(CI) 



where <f „ = (^n)ai3js(p, q, P) is the Dirac doublet in 
Eq. (Bl) involving the mixed-antisymmetric and 

mixed-symmetric (^2) contributions, and the rep- 
resent the three diagrams in the Faddeev equation (2). 
We write the amplitude decomposition as 



*„(p,</,P) = X;/„,<(t)X i (r,«,P), 



(C2) 



where the index i now collects all the previous indices 
{k,i,j,(jj} from Eq. (BIO), and t = {p 2 , q 2 , Zq, Z\, z-i\ 
again abbreviates the invariant momentum variables. Us- 
ing the trace orthogonality (B16) yields equations for the 
amplitude dressing functions: 



/»,<(*) = £/#(*), 

a=l 

where the contributions on the r.h.s. are given by 



(C3) 



(C4) 



Here the g£\ denote the coefficients of the wave functions 

$^ a ) = S{pb)S(p c )'^ n with attached propagator legs, 
where {a, b, c} is a symmetric permutation of {1,2,3}, 
with the same Poincare-covariant decomposition as in 
Eq. (C2). The internal momentum variables v- a > are ob- 
tained from the internal momenta of Eq. (5). The kernel 
and propagator matrices in (C4) follow from a projection 
onto the basis elements: 



X K aa >ppi (k) , 

9i? = I *4 (r, s, P)p a s~f Xj (r, s, P) a >p> 7 6 



(C5) 



X S aa >(pi) Spi3'{p2) 



and a cyclic permutation of the index pairs and quark mo- 
menta in the second and fourth row yields the remaining 
expressions for a = 1, 2. 

Eqs. (C3-C4) can be solved by iteration which requires 
a compromise between huge memory requirements for the 
kernel and propagator matrices (if computed in advance) 
or large CPU times (if computed in each iteration step 
anew). Compared to the permuted diagrams for a = 1,2, 
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the kernel K,^ is less difficult to handle: in the kinemat- 
ics of Eq. (5) the internal relative momentum p^> equals 
the external one p, and hence A 3 ^ = r which in the rest 
frame reduces to e%. In combination with the simplicity 
of the rainbow-ladder kernel which only depends on the 
gluon momentum k — — q, the kernel K.^ is indepen- 
dent of the variables p 2 and Z\ and hence only requires a 
modest amount of memory (typically several GB). 

Nevertheless the storage problem for K.^ and K.^ re- 
mains. In this respect it is beneficial to take advantage 
of the permutation-group properties of the amplitude. 
Applying the relations of Eq. (B5) for the internal am- 
plitudes #(pW,gW,P) and f(?/ 2 \ q {2 \ P) in the Fad- 
deev equation (2), and expressing the relative momenta 
in each diagram by p' , q' or p" , q", respectively, yields 
the relation 



q, = #(3) + M ' [* (3) ] ' + M" [* (3) ] ' 



(C6) 



By means of Eq. (C6), the amplitudes and no 
longer need to be computed explicitly but can be recon- 
structed from \E'( 3 ). Note that the full Dirac-color- flavor 
amplitude corresponding to is not totally antisym- 
metric, i.e. its Dirac parts do not transform via Eq. (B5), 
and the right-hand side of the above equation is not sim- 
ply 3*^ . 

The orthogonality relation (B16) then yields corre- 
sponding equations for the amplitude dressing functions: 

/<(*) = /i 3) (*) + M'HL jf (0 + M"H!L /, (3) (t") (C7) 



where t' and t" are the Lorentz-invariant momentum vari- 
ables obtained from the permuted momenta (B4); M! 
and M." act on the doublet index n; and the matrices H' 
and H" are given by 



; j (t) = iTr{x j (r, S ,P)X;.(r', S ',F)}, 
l ^t) = \Tr{% i {r,s,P)Y!l{r",s",P)) 



(C8) 



where Xj , X" are the basis elements with permuted Dirac 
indices according to App. B 1, and the unit vectors r' , r" , 
s' , s" follow from the permuted momenta in Eq. (B4). 
Upon exchanging p 2 and q 2 with the singlet and doublet 
variables x and y\ of Eq. (B7) one finds that H' and 
H" arc independent of x; and by taking into account the 
sparseness of )C\^ and H', H" , these matrices can be 
stored in advance. Instead of interpolating the dressing 
functions for the permuted momentum variables in each 
iteration step again, one may expand them in orthogonal 
polynomials. Due to the weak angular dependencies in 
zq, Z\ and z^ (and also yi) sufficient accuracy is already 
reached by retaining a small number of moments. 

The algorithm for solving the Faddeev equation can 
now be summarized as follows: start with a guess for 
fn,i{t)l m each iteration step, compute /„/(£) from 
Eq. (C4) and f n ,i(t) from Eq. (C7); proceed until con- 
verged. These steps are repeated for different nucleon 



test masses for which purpose an eigenvalue A(M ) is in- 
troduced in the equation; the correct nucleon mass yields 
the eigenvalue A(M) = 1. With the procedure detailed 
here, solving the Faddeev equation on a cluster (or even 
a potent desktop computer), without any approximation 
on the momentum dependence and with accuracy that is 
sufficient to determine a reliable nucleon mass, becomes 
possible within a few hours. 



Appendix D: Nucleon electromagnetic current 

1. General properties 

The matrix-valued electromagnetic current of the nu- 
cleon can be written in the most general form as 



"2) iT 



P 2 



~M 



AL 



(Dl) 



where /i is the photon index and Jg, s (P,Q) is a Dirac 
matrix with indices S, 6' for incoming and outgoing nu- 
cleon amplitudes. We abbreviated the positive-energy 
projectors by A + (Pf^) — A^'\ The current involves 
two momenta, expressed through the incoming and out- 
going momenta Pj, Pt or by the average momentum 
P = (Pi + Pf)/2 and photon momentum Q = Pf — Pi. 
Since the nucleon is on-shell, P 2 = P 2 = —M 2 , one has 

P 2 = -M 2 (l + t), P-Q = 0, (D2) 

where r := Q 2 /(4Af 2 ); hence the Lorentz-invariant 
form factors which constitute the vertex can only de- 
pend on the photon momentum-transfer Q 2 . Contracting 
Eq. (Dl) with nucleon spinors u(P/,s/), u(P,-,Sj) which 
are eigenstates of A+ l , e.g.: 



A t ,u(P i ,s) = u(Pi,s) , 



±1/2. 



(D3) 



yields the current-matrix element (Pf, Sf \ J M | Pi, Si). 

As a three-point function depending on two momenta, 
with one vector and two spinor legs, the bracket in 
Eq. (Dl) could in principle involve the 12 positive-parity 
tensor structures 

{r,P»,Q»}x{t,f, Q, [f,Q]}. (D4) 

Applying the positive-energy projectors reduces this set 
via Eq. (All) to the three basis elements { 7 M , P M , Q M }, 
and imposing charge-conjugation invariance of the cur- 
rent, Jv(P,Q) = J fi (P,—Q), eliminates the component 
Q M . The resulting current is automatically conserved, i.e. 
Q^Jf* = 0. Consequentially, the general electromagnetic 
current of the nucleon depends only on two form factors 
Fi(Q 2 ). Using the Gordon identity 

iP» a^ v Q v ' 
~M + 2M 



r 



M = 



(D5) 



finally leads to the expression given in Eq. (18), and the 
electric and magnetic form factors are extracted via the 
Dirac traces in Eq. (33). 
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2. Diagrams in the three-quark framework 

In this appendix we collect the ingredients of the nu- 
cleon electromagnetic current operator in the three-quark 
framework which is depicted in Fig. 7. In a rainbow- 
ladder truncation the current is the sum of the impulse- 
approximation and kernel diagram, the sum of the three 
permutations a = 1,2,3, and the sum of the mixed- 
antisymmetric and mixed-symmetric components in the 
incoming (n = 1,2) and outgoing (n 1 = 1,2) nucleon 
amplitudes: 



^(^)=EE[^ IMp +e 



a—l n'n 



I-' 

S'S 



(D6) 



For example, the impulse-approximation diagram where 
the photon couples to the upper quark leg reads explic- 
itly: 



{y n ')p> a >8>y(Pf,qf,Pf) x 

(D7) 
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where pi, qi and pf, qf are the incoming and outgoing 
relative momenta; Pi and Pf are the incoming and out- 
going nucleon momenta; pi, P2, P3 andpf — p 3 ±Q/2 are 
the quark momenta; S(pt) are the dressed-quark propa- 
gators; is the dressed quark-photon vertex; and p and 
q are the two loop momenta. For symmetric momen- 
tum partitioning (i.e., a momentum-partitioning param- 
eter 1/3) the relative momenta are explicitly given by 



Pf 



, Q 
P + J 



Pi=P~ V 



Q 
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and the quark momenta by 
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and p 3 = p + P/3. 

The flavor trace in Eq. (D7) denotes 
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where Q = diag(q u , qa) is the quark charge matrix at- 
tached to the quark-photon vertex and F n , FL are the 
flavor matrices of Eq. (B2). Applying the isospin vec- 
tors ei = (1,0) or &2 — (0,1) in Eq. (D10) singles out 
the proton and neutron contributions, respectively. This 
yields 



p(3) _ f qu 

for the proton and 
F (3) = 







\ (q u + 2q d ) 
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for the neutron flavor traces. The color traces for the 
impulse- approximation diagrams equal —1. 

The current diagram involving the rainbow-ladder ker- 
nel has the same overall shape as Eq. (D7), 
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x Sa-MSppfa) [S(pi)T^p 3 ,Q)S(p^)}^x 

X (*l 3) )a/3 T 5(Pi,%,-Pi), 
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except for a color factor 2/3 and an incoming amplitude 
which is replaced by: 

(^)<xp-ys(Pi,qi,Pi) = / K aa ,p p ,(k) x 



X S a > a "(Pl ~ k) Sp'i3"{p2 + k) (^n)a"0"j6{Pi,qi + k,Pi) . 

The practical form-factor calculation can be simplified by 
writing the sum of Eqs. (D7) and (D13) schematically as 
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where <Fj is the wave function SiS 2 S 3 ^i, and $| 3 ' ) the 
wave function obtained from the third diagram in the 
Faddeev equation. Both of them can be collected before- 
hand when solving the equation and are implemented in 
the form factor diagram simply through evaluation at the 
proper incoming momenta, such that only the product 
53 r^ 1 (which is independent of the loop momentum q) 
needs to be computed explicitly. 

Again the diagrams for a = 1,2 can be inferred from 
the a = 3 diagram through permutations, namely: 
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To prove this, one writes down the expressions for 
j(o=l,2) ana i g OUS to Eqs. (D7) and (D13) with respec- 
tive momentum dependencies on p, q, P and Q; exploits 
the doublet transformation properties (B5) for the am- 
plitudes; expresses all internal momenta through p' and 
q' (for a = 1) or p" and q" (a — 2) from Eq. (B4); and 



replaces the 1 



integration over {p, q} by integrating 



over {p',q'} or {p",q"}, respectively. 

Evaluation of Eq. (D14) finally yields the following ex- 
pressions for the total nucleon electromagnetic current of 
Eq. (D6): 
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If we denote by G AA the form factors obtained via 
Eq. (33) from the matrix element J{ x involving the 
mixed-antisymmetric Dirac amplitudes, and by G 
those corresponding to the mixed-symmetric overlap 
J 22 , one obtains at the relation (35) for the proton and 
neutron form factors: 



G p = 2G 
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G n = G 



ss 
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To conclude this section we note that we encountered 
difficulties in the form factor calculation in connection 
with the angular variable z 2 of Eq. (12). The dependence 
upon z 2i = % ■ Pi and #2/ = • Pt in the incoming 
and outgoing nucleon amplitudes of Eqs. (D7) and (D13) 
is reconstructed from a Chebyshev expansion in the rest 
frame as obtained from the Faddeev equation. Such a re- 
construction is sensitive to the numerical accuracy in the 
Faddeev amplitude and, especially for z 2 , leads to con- 
vergence problems at larger Q 2 . To some extent this is 
expected, as \z 2 i\, \z 2 /\ < \f\+~r in the form-factor inte- 
gral, where r = Q 2 / (4M 2 ), i.e. the domain of the expan- 
sion is no longer bounded by a unit circle and grows with 
the photon momentum. Eventually, with more powerful 
computing resources, such issues can be avoided alto- 
gether by solving the Faddeev equation in each moving 
frame anew, cf. Ref. [60] in the context of the pion. While 
the dependence on z 2 in the rest frame is very weak and 
could be neglected in the form factor calculation, its in- 
clusion is still necessary to avoid issues with charge con- 
servation at Q 2 — 0. To account for this we perform the 
full Chebyshev resummation in z 2 i, z 2 f at each Q 2 only 
if the conditions 



Pi. 
4 



> Q 2 , x f :-- 
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Pf q f ^2 
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are satisfied. By virtue of Eq. (D17), the full angular 
dependence is taken into account at Q 2 — 0, whereas at 
large Q 2 the photon momentum acts as a cutoff on the 
Chebyshev expansion and only the zeroth moments are 
retained. This procedure resolves the problem mentioned 
above. 



3. Quark-photon vertex 

The only ingredient in the form factor diagrams of 
Fig. 7 which has not already been defined in connection 
with the Faddeev equation is the dressed quark-photon 
vertex T^(k, Q). Its general expression is derived from 
the Ward-Takahashi identity 



Q^mk,Q) = S- 1 (k + )-S- 1 (k_) 
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and by imposing regularity at Q 2 = 0. It is expressed by 
a sum of the Ball-Chiu term [118] and a purely transverse 
contribution: 

r£(fc, Q) = ^ A + ^iik 1 A A + A s ) + T£ , (D19) 
where k± = k ± Q/2 are the incoming and outgoing 
quark momenta and A(k 2 ) and B(k 2 ) — M(k 2 )A(k 2 ) are 



the dressing functions of the inverse quark propagator 
S-^k) = i^A(k 2 ) + B(k 2 ), with 



S P := 



F(fci) + F{k 2 _) 



A P := 



F(k 2 
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F(k 2 _] 



The transverse part can be written as 
-iT» T = 1 »(h+if 2 ®) + 

+ i/a k-Q\ [<#, n + U\ b£, fa] Q + (D20) 
+ k£(if 5 + f 6 k-Q@ + f 7 k' + if s fa Q) , 

where the fi(k 2 , k ■ Q, Q 2 ) are scalar dressing functions 
and 7^, fc^, are transverse with respect to the photon 
momentum Q. While this basis decomposition is not 
orthogonal, it provides a simple representation for the 
Ball-Chiu vertex if it is used to describe the vertex in 
total. 

The quark-photon vertex is obtained self-consistently 
from its inhomogeneous Bethe-Salpeter equation [114], 
given by 



a/3 



4 f n (D21) 

Krtpfi [S(k' + )T^k',Q)S(k'_)] a , , , 



where Z 2 is the quark renormalization constant and 
Kaa'pip the rainbow- ladder kernel of Eq. (8) that also 
appears in the quark DSE and the covariant Faddeev 
equation. 

To solve the inhomogeneous BSE numerically it is con- 
venient to express the vertex in an orthonormal basis, for 
instance 

r?(r,Q) G {^7f, r", x {l, 0, f, /#} (D22) 

where, in analogy to the nucleon amplitude decompo- 
sition of Eq. (BIO), we defined a relative momentum 
transverse to the photon momentum by the unit vector 
r := fey, and in addition a 7— matrix which is transverse 
to both Q and r, i.e. 7^ := 7^ — r^f. The basis elements 



satisfy the orthogonality relation 



(D23) 



Upon projecting the kernel onto this basis it becomes 
clear that the equations for the transverse (~ 7f,r M ) 
and longitudinal parts (~ Q M ) decouple and can be 
solved independently. The longitudinal result reproduces 
the longitudinal projection of the Ball-Chiu vertex and 
thereby the Ward-Takahashi identity, and purely lon- 
gitudinal terms do not contribute to nucleon form fac- 
tors because of current conservation: Q M J M = 0. Hence 
it is sufficient to consider the eight transverse elements 
alone. The inhomogeneous BSE self-consistently gener- 
ates a timelike vector-meson pole in the quark-photon 



vertex at Q = 
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